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Abstract 

We demonstrate that human electrophysiological recordings of the local field potential (LFP) from in- 
tracranial electrodes, acquired from a variety of cerebral regions, show a ubiquitous 1/f 2 scaling within the 
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power spectrum. We develop a quantitative model that treats the generation of these fields in an analogous 
way to that of electronic shot noise, and use this model to specifically address the cause of this 1/f 2 Brow- 
nian noise. The model gives way to two analytically tractable solutions, both displaying Brownian noise: 1) 



>• uncorrected cells that display sharp initial activity, whose extracellular fields slowly decay and 2) rapidly 
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ON firing, temporally correlated cells that generate UP-DOWN states. 
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I. INTRODUCTION 

Power laws appear in a large variety of settings throughout nature and often signify that there 
is a simple process at the origin of what appears to be a very complex phenomena. Examples of 
the variety of settings in which power laws appear are the Gutenberg-Richter law for the size of 
earthquakes [1], the allometric scaling laws that appear throughout biology 0, and Paretos's law 
of income distributions J3]|. 

Power laws have also been witnessed within the brain in electroencephalographic (EEG) and 
magnetoencephalographic (MEG) recordings while studying a wide variety of brain function flU. 
The signals recorded outside the skull by these techniques represent the global activity of a large 
amount of cortical and subcortical tissue and give rise to a range of exponents. Much more local 
measurements of cerebral activity may be recorded by a single microelectrode. While the em- 
phasis of these measurements is usually focused on the spiking activity of single cells within the 
vicinity of the recording electrode, local field potentials (LFPs), which comprise the much slower, 
background of electrical activity, may also be extracted from the signal. 

We found that electrophysiological recordings, taken from pharmacologically intractable 
epilepsy patients with microelectrodes implanted in a variety of cerebral areas, display a sur- 
prisingly universal l// 2 power law in the frequency spectrum of LFP activity (Fig. [I]). A power 
spectrum of this sort is said to display the statistics of Brownian noise since it has the same scaling 
exponent as a ID random walk. However, it is far from clear what the underlying mechanism is 
that gives rise to these statistics. 

To address this question, we developed a general method for modeling the LFP from what we 
refer to as neuronal shot noise, that allows one to study the microscopic origin (i.e., single neuron 
activity) of the power law dependence in the power spectrum. We propose two quite different 
processes that could both give rise to the observed 1/ f 2 dependence. The first involves the uncor- 
rected firing of single neurons that display very slow dendro- synaptic decay in the extracellular 
field which they generate. The second possibility involves the correlated firing of a single neuron 
which displays either no activity (DOWN state) or very rapid spiking (UP state). We end with 
a discussion of the UP-DOWN states (UDS) suggested by our model and how they compare to 
experimentally observed UDS within the cortex. 
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FIG. 1 : Power law and scaling exponent in local field potentials recorded from the human cerebral cortex. 
A. Exemplary power spectrum of local field potentials recorded from a micro-wire in the temporal lobe. 
The scaling exponent (here a = 2.04) was determined by a linear least-square fit of the log-log power 
spectrum (see inset). B. Scaling exponents (mean ± stand, dev.), averaged across micro-wires for different 
brain regions. FL: frontal lobe; TL: temporal lobe; PL: parietal lobe; OL: occipital lobe. 

H. BROWNIAN NOISE IN HUMAN ELECTROPHYSIOLOGY 



We recorded local field potentials from the cerebral cortex of 10 subjects with pharmacolog- 
ically intractable epilepsy (4 males; 24 — 46 years old), implanted with chronic electrodes to 
localize the seizure focus for possible surgical resection. Electrode locations were based exclu- 
sively on clinical criteria and were verified by MRI or by computed tomography co-registered to 
preoperative MRI. Each electrode probe had nine micro-wires (Platinum/Iridium, 40/im diameter) 
protruding from its tip, eight high-impedance recording channels (typically 200 — 400kfi) and one 
low-impedance reference with stripped insulation. The differential signal from the micro-wires 
was amplified using a 64-channel Neuralynx system, filtered between 1 and 9000 Hz, and sam- 



pled at 28 kHz. We recorded from a total of 684 micro-wires (106 in the frontal lobe, 546 in the 
temporal lobe, 16 in the parietal lobe, 16 in the occipital lobe). Recordings lasted for 10 min while 
subjects were awake and at rest with eyes open. All studies conformed to the guidelines of the 
Medical Institutional Review Board at UCLA J5]|. 

For analysis, the data was down-sampled to 7 kHz using an anti-aliasing filter. The power 
spectral density was estimated by applying Welch's method to consecutive 5-sec segments and 
subsequently averaging over the entire 10 min (Fig. 1A). The scaling parameter a was determined 
as the slope of a least-square linear fit of the double-logarithmic power spectrum. To diminish the 
influence of amplifier roll-off, the linear fit was restricted to a frequency range of 1 to 400 Hz (Fig. 
1A, inset). Figure IB shows the scaling parameters averaged across different micro-wires for four 
different regions of the cerebral cortex along with their standard deviation. Note that in all four 
regions the scaling parameter is close to a value of a = 2, indicating a universal scaling behavior 
of local electrical brain activity. To the best of our knowledge, this is the first time such a universal 
feature of the LFP has been reported in humans. 

III. NEURONAL SHOT NOISE 

The microscopic generation of the Local Field Potential (LFP) may be formulated in a similar 
way to that of shot noise, originally described by Schottky J6l to account for the noise across an 
electrical resistor. This may be seen by writing the extracellular potential V(t) generated by TV 
neurons, at a given spatial location within the brain, as follows: 

N 

V(t) = J2 dt'fiit-t')^), (1) 

i J 

which is exactly how one quantitatively models shot noise. Here the function fi(t) accounts for the 
temporal profile of the extracellular field generated by neuron i while \Xi (t) represents the activity 
of that neuron and may be explicitly written as Hi(t) = J2k $(t — t k .), where 5(t) is the Dirac delta 
function. From this definition, we see that the function /Xj(t) may be thought of as analogous to 
the spike train with firing activity occurring at times t^ for neurons i = 1 . . . N. Note, this model 
does not require that the neurons generate action potentials, it only assumes a stereotyped profile 
fi(t) for the electrical field generated by each neuron which repeats at times t^ (see Fig. [2]). 

We will assume that the relevant neural activity has reached a steady- state such that the auto- 
correlations G(t) = (V(t)V(t + r)) are independent of t. By the Wiener-Kinchin theorem, the 
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FIG. 2: Schematic of Eq.fTlrepresenting neuronal shot noise. A function fi(t), representing the extracellular 
field associated with the i th neuron, occurs at times t\, ti . . . governed by the statistics of Hi(t). 

autocorrelation function G(t) is simply related to the power spectrum S(cu) by Fourier transfor- 
mation. From Eq. [T] we may write the power spectrum as 

S(u) = (\V(uj)\ 2 ) ^ftWfjMfaMiijM), (2) 

hi 

where V(u>), f(ui), and jx(u) are the Fourier transforms of their respective temporal functions. 



IV. BIOPHYSICAL EXAMPLES 

To solve for the power spectrum as written in Eq. [2] would require us to know the location of 
each neuron involved in generating the LFP, the extracellular field produced by each neuron, and 
the decay of that field through the neuronal medium. While we have carried out such biophysical 
calculations in the past for single neurons 0, we are here only concerned with understanding the 
source of the 1/ f 2 behavior of the power spectrum, not in reconstructing the LFP. 

The scale invariant nature of S(u) greatly simplifies our problem since we need not concern 
ourselves with any constants that arise from the details mentioned above. To clarify this point, 
let us assume that the power spectrum S(cu) = C\uj n , where C\ is a constant. We can solve 
for the coefficient of the power by plotting the log of both sides of this equation, log S(uj) = 
n logo; + log C\. The power dependence is given by the slope n and is unaffected by the constant 
offset. 

Since the power law dependence of -S'(a;) is independent of overall constants, we may neglect 



the previously mentioned biophysical details and focus on finding solutions to Eq. [2] that have a 
1/uj 2 functional dependence. In general, this problem is still quite difficult; however, there are two 
limits which allow a simple, analytical solution. We will now discuss these two cases. 

Case I: Slow Dendro-Synaptic Decay. The simplest case to consider is that the spiking statistics 
are independent between neurons, and that the spiking of each neuron is an independent Poisson 
process |[8|. In this case 

(jii^Hjiw)) = 6 id p$, (3) 

where fy is the average firing rate of the neuron i and 5 id is the Kronecker delta function. 
We can now ask when the power spectrum satisfies 

sm = 5> 2 i;,mi 2 oc4 (4) 

i 

The solution requires fi(u) oc 1/uj, whose Fourier transform is a Heaviside step function fi{t) oc 
9(t). This answer is a bit unrealistic since it implies that the field generated by the cell does not 
decay. A more realistic solution would be to assume a form such as 

fi(t) oc 6(t)e- a \ (5) 

which has Fourier transform I /(a + iuj). In the limit of slow decay, a< 1, a neuron with an 
extracellular field of this form, firing with Poisson statistics, would give rise to Brownian noise in 
the LFP. 

In this case, the l// 2 behavior originates from the steep onset of the extracellular field. The 
rise time of an action potential may occur within a fraction of a millisecond, which could account 
for a sharp onset, while the decay of the dendro- synaptic extracellular field may last for as long 
as a second [9J. The functional form of the decay does not affect these results, so long as the cell 
takes much longer to return to baseline than it took to spike. 

Case II: UP-DOWN States. The second case that we consider is the limit of a sharply peaked 
extracellular field. In this case, we may treat fi(u) as a constant / is and we will assume that the 
activity of different neurons are either uncorrelated ((/ij/ij) oc 8 id ) or synchronous ((/Xj/Xj) oc 1). 
The spike timing of a single neuron, however, may show a temporal correlational structure. These 
assumptions lead to a power spectrum 
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FIG. 3: Top: Binary sequence generated from a telegraph process. Bottom: Power spectrum of the binary 
sequence confirming the l// 2 behavior. The slope of the dashed line is -2. 

Since all the frequency behavior is contained within the statistics of fa, and we are assuming that 
all cells are active with the same statistics, we need to look for a sequence of spikes that have 
individual spike timing correlations of the form 



(//(£ + r)/i(t)) oc r. 



(7) 



since linear time correlations are consistent with l/cu 2 frequency correlations. This is the same 
linear in r scaling as that of a ID random walk and is at the origin of the term Brownian noise. 

Since /j,(t) is analogous to the spike train of each neuron, we need to formulate a binary se- 
quence that shares the correlational structure of a random walk. A simple way to generate a binary 
sequence representing white noise is to pick a random number at each timestep and then apply a 
threshold such that all values above the threshold are set to one, and all below to zero. Brownian 
noise may be created by integrating a white noise signal. However, it is not obvious how to apply 
a similar thresholding procedure to convert the resulting analog signal into a digital one [10]. 

An alternative procedure that will generate Brownian noise is given by setting up a telegraph 
process [11J. In this case a binary sequence is generated by constructing a two-state system, (0 



and 1) where, at each timestep, the probability of making the transition — » 1 (1 — ► 0) is given by 
k + At (/c_ At). The autocorrelation function for such a process may be explicitly written as 

M^M*)H(^) 2 (i + ^-<^). (8) 

In the limit of equal transition rates, k + = k—, and low probability of making a transition, [k + + 
A;_)r <C 1, Eq. [8]reduces to 

(M* + tM*))«^(1-*t), (9) 

which has the desired linear in r statistics of a random walk. 

Figure [3] displays a binary sequence generated by a telegraph process and the 1/f 2 dependence 
of its power spectrum. The telegraph process gives rise to periods of sustained, rapid activity 
followed by intervals of inactivity. This results in collective oscillations that display a much lower 
frequency than the rapid firing witnessed during depolarization . The result is a pattern of behavior 
reminiscent of UP-DOWN states common in cortex lfT2l . 

V. DISCUSSION 

Studies of the LFP and single neuron spiking activity, combined with current source-density 
analysis, suggest that LFPs are primarily the result of dendritic activity distributed over a large 
region of the cortex. LFPs are therefore believed to provide a measure of the local processing and 
input to a given region of the brain [fT3l . 

We developed a very simple model to explain our key experimental finding, a 1/f 2 decay in the 
local field potential recorded from individual microelectrodes implanted into human cortex. In par- 
ticular, we showed two examples of how biologically realistic, single neuron activity, parametrized 
by the temporal shape of the resulting extracellular fields and the statistics of cellular activity, can 
give rise to power law behavior within the LFP. 

In Case I , we showed how a population of cells, each displaying a sharp onset of activity and 
a much slower decay of the extracellular field, could give rise to a Brownian power spectrum. 
The time course of dendritic activity is often much longer than that of an action potential. This 
is in line with the above statement concerning the origin of the LFP. However, the sharpness of 
the temporal onset of activity is what gives rise to the power-law behavior. One mechanism that 
might account for this result would be the rapid initiation of an action potential, followed by slow 
dendro-synaptic decay. 

8 



Of course, this model not only assumes that the spiking statistics of each neuron is Poisson, 
but that there are no correlations among neurons. It is not uncommon to find the firing rate of 
single neurons uncorrected with the averaged behavior of the local population, however, this is 
not always the case |fl4l. 

In Case II, we explored the opposite extreme from Case I, that of rapidly firing, single neurons 
with linear temporal correlations. This behavior is similar in nature to so-called UP-DOWN states 
seen in cortical neurons. During periods of sleep, quiet awake behavior, or under a variety of 
anesthetics, spontaneous activity of neocortical neurons display slow 0.1 to 2 Hz oscillations called 
UP-DOWN states (UDS). These states appear to be characteristic of slow-wave sleep 03] an d are 
thought to be involved in the consolidation of long-term memories and in learning. The UDS of 
cortical pyramidal neurons are highly synchronized and may be clearly seen in LFP recordings 
of the cortex. The UP states are characterized by a sustained depolarization that leads to rapid, 
20-70 Hz spiking activity while the DOWN states show periods of hyper-polarized inactivity. 

It should be pointed out that our recordings were performed in the awake resting state in the 
human cortex, whereas UDS and ultra-slow oscillations have been described only in states of low 
vigilance such as slow- wave sleep and anesthesia in animal studies. It is therefore unlikely that 
the power law scaling behavior observed in our recordings would be caused exclusively by the 
mechanisms illustrated in Case II. Nevertheless, it is encouraging that this extreme analytical case 
of a 1/f 2 power law scaling gives rise to phenomena that are actually observed in mammalian 
brains. 

The true origin of the l// 2 behavior probably lies somewhere in-between the two limiting 
cases we have considered here. Unfortunately, an analytic evaluation of Eq. [T] when there is an 
explicit time dependence in both the extracellular field (fi(t)) and the firing statistics (/ij(t)) is, in 
general, difficult. However, for a known set of fi(t) and /ij(t), a numerical evaluation of Eq.[T]is 
straightforward. This formalism should, therefore, serve as a starting point in modeling power-law 
dependencies in the power spectrum of the LFP and in connecting this property to the underlying 
single neuron activity. 
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